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ABSTRACT 

Context. The cosmic infrared background (CIB) consists of emission from distant, dusty, star-forming galaxies. Energetically, the 
CIB is very important as it contains as much energy as the extragalactic optical background. The nature and evolutionary status of the 
objects making up the background are, however, unclear. 

Aims. The CIB peaks at ~150yum, and as such is most effectively studied from space. The limited apertures of space-borne tele- 
scopes set the angular resolution that can be attained, and so even Herschel, with its 3.5m diameter, will be confusion-limited at this 
wavelengths at ~5mJy. The bulk of the galaxies contributing to the CIB are fainter than this, so it is difficult to study them without 
interferometry. Here we present the results of a preliminary study of an alternative way of probing fainter than the continuum confu- 
sion limit using far-IR imaging spectroscopy. An instrument capable of such observations is being planned for SPICA - a proposed 
Japanese mission with an aperture equivalent to that of Herschel and more than 2 orders of magnitude more sensitive. 
Methods. In this paper we investigate the potential of imaging spectrometers to break the continuum confusion limit. We have 
simulated the capabilities of a spectrometer with modest field of view (2'x2'), moderate spectral resolution (R~ 1-2000) and high 
sensitivity. 

Results. We find that such an instrument is capable of not only detecting line emission from sources with continuum fluxes substan- 
tially below the confusion limit, but also of determining their redshifts and, where multiple lines are detected, some emission line 
diagnostics. 

Conclusions. 3-D imaging spectrometers on cooled far-IR space telescopes will be powerful new tools for extragalactic far-IR as- 
tronomy. 
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1. Introduction 

The discovery of the Cosmic Infrared Background (CIB) (Puget 
et al., 1996; Fixsen et al, 1998) provided a key step in our un- 
derstanding of the history of the universe. The energy in the 
CIB is comparable to that in the integrated UV/optical/near-IR 
(UVOIR) background. Since the CIB comes from obscured sys- 
tems and the UVOIR from unobscured systems this means that 
roughly 50% of energy produced in the history of the Universe 
comes from obscured systems. The nature of the objects mak- 
ing up the CIB, however, remains unclear. Source populations 
that contribute significantly to the CIB have been found in sub- 
millimetre surveys (Blain et al., 2002, and references therein), 
in the far-IR (Puget et al., 1999), and in the mid-IR (Dole et al., 
2006), with as much as 70% (Dole et al., 2006) of the CIB now 
believed to have been resolved into individual sources. Much of 
this work, however, has been done at wavelengths far away from 
the 150yt/m peak of the CIB, and is reliant on extrapolations. 
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The Herschel far-IR mission (Pilbratt, 2003), to be launched 
in 2008, will be capable of deep observations around the peak 
of the CIB, and will be sensitive enough to resolve a signifi- 
cant fraction of CIB sources at wavelengths where their spec- 
tral energy distributions (SEDs) peak (Dole et al., 2004). With 
a 3.5m-diameter mirror, Herschel will have the largest satellite 
primary to date, but even this will only give an angular resolu- 
tion of ~10" at 150/im. The source density of CIB sources on the 
sky is such that confusion, rather than instrument sensitivity, is 
the main limitation on the faintness at which an individual CIB 
source can be detected - for exmple, the confusion limit of the 
PACS instrument on Herschel (Poglitsch et al., 2005) has been 
estimated to be ~5mJy at the ~ 150jt/m peak of the CIB (Jeong et 
al., 2006). Extrapolations from the deep SPITZER surveys sug- 
gest that Herschel will therefore only be able to resolve ~50% of 
the CIB (Dole et al., 2004) at best, even at its shortest operating 
wavelength. 

There are of course ways in which one can probe below the 
confusion limit by using observations at other wavelengths or 
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through statistical approaches. This includes extraction of long 
wavelength fluxes at the positions of sources detected at shorter 
wavelengths which has proved successful in some studies (eg. 
Dole et al., 2006). There remains the worry that such studies 
can be biased against redder sources, which may be missing 
from the short wavelength catalogs but still contribute strongly 
at long wavelengths. Such issues can only be avoided by detect- 
ing sources at the long, CIB, wavelengths themselves which will 
require observations below the continuum confusion. In this re- 
search note we put forward the results of work that has been 
done to use an additional observational dimension to break the 
confusion limit - that of spectroscopy. We will demonstrate that 
through spectroscopy one can identify sources that lie signifi- 
cantly below the traditional (contiuum) confusion limit, and thus 
to probe a population that contributes significantly to the CIB, 
that cannot presently be studied in any other way. 
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Table 1. Galaxies used to provide far-IR lines for templates. 



2. Modeling Blank Field Far-IR Line Surveys 

2.1. A Far-IR Imaging Spectrometer 

The key feature of an instrument designed to conduct blank field 
searches for emission lines is that it must simultaneously cover 
a significant field of view over a significant spectral range with 
sufficient sensitivity to detect the lines against the continuum. 
Imaging spectroscopy at far-infrared, infrared and optical wave- 
lengths can be achieved using a number of different instrumental 
techniques, including imaging fourier transform spectroscopy - 
as employed in the Herschel SPIRE instrument (Griffin et al., 
2004), integral field spectroscopy using a grating and image 
slicer - as employed in the Herschel PACS instrument (Poglitsch 
et al., 2005), using imaging Fabry-Perots - as employed in SPIFI 
(Stacey et al. 2002), and at the longer wavelengths with het- 
erodyne arrays using ultra-sensitive mixers - as employed in 
HARP-B (Smith et al., 2003). Each technique is more or less ap- 
propriate to a particular application and a particular wavelength 
range/instantaneous bandwidth requirement. 

Instrument concepts based on the first two of these tech- 
niques are currently under consideration for an imaging spec- 
trometer (ESI, the European SPICA instrument) for the 35- 
210,1*111 waveband for SPICA (Swinyard et al., 2006). For the 
purpose of this paper our considerations of the instrument it- 
self are restricted to its deliverables, namely a data cube with a 
particular two-dimensional spatial extent and angular resolution 
and a given spectral range and resolution in the third (spectral) 
direction, with a specified sensitivity in each spectral channel. 
We have based our simulations on what we believe to be feasi- 
ble for a cooled ( 4K), 3. 5m-diameter primary, telescope such as 
SPICA: a spatial resolution of 8" as set by the diffraction limit 
at ~120//m; an instantaneous spectral coverage of 60 to 210/im 
and a fixed channel width of AA * 0.176yum, resulting in a spec- 
tral resolution ranging from 350 at 60yum to 1200 at 210^111 and 
an instantaneous field of view of 128"xl28". These instrumental 
parameters are very similar to those expected in the two longest 
wavelength channels of the ESI far-infrared instrument concept. 

2.2. The Far-IR Sky 

To assess the ability of our far-IR imaging spectrometer to de- 
tect line-emitting sources in the CIB, we must model the far- 
IR sky. Our starting place is the current generation of empirical 
number count models that fit all current constraints on galaxy 
counts from optical to submm, including the strength of the CIB. 
The model we have chosen to adopt is the burst mode evolu- 



tion model of Pearson (2006) since this provides all the neces- 
sary features we require, and, in later work, can easily be com- 
pared and contrasted with Pearson's parallel bright galaxy evo- 
lution model. The model we use is based on an evolving, type- 
dependent, far-infrared luminosity function which makes use of 
a set of 1 1 optical-to-submm template SEDs. The template ob- 
jects range from normal, Milky- Way type galaxies, through in- 
creasingly far-IR luminous objects to extreme ultraluminous in- 
frared galaxies (ULIRGs) and type 1 and 2 AGNs. To simulate a 
region of the far-IR sky we begin with a simulated catalogue of 
sources with far-infrared fluxes, luminosities, redshifts and real- 
istic SEDs taken from the galaxy evolution model. Each source 
in the catalogue is parameterised by a wavelength and a flux, 
which are assigned based on the SED template type. At this 
point, the template fluxes are continuum only. Far-IR emission 
(and absoprtion) lines are added to these templates using line 
strengths derived from LWS far-IR spectroscopy for a range of 
nearby galaxies (eg. Negishi et al. (2001)), as tabulated in Table 
1 . The resulting templates, including lines, are shown in Figure 
1 . We only include the brightest far-IR lines, and no MIR lines, 
at this stage since these will be the most useful for source detec- 
tion in the far-IR. 

To generate the data cube that will come from observing this 
catalog we compile a list of sources in our field of view. The inte- 
grated flux contributed to each 8" pixel by all the objects that lie 
within it is then calculated. The template spectrum for each ob- 
ject contributing to a spatial pixel is shifted in wavelength by the 
source redshift and then normalised to the flux of the source at 
the fiducial catalog wavelength (40//m, the wavelength at which 
the original catalog was normalised). The flux at a given, ob- 
served wavelength is then interpolated from the template and 
added to the appropriate wavelength channel in the data cube. 
This process is repeated for each 8"x8" spatial pixel to build up 
the 16 x 16 spatial cube, each pixel made up of a 3rd dimension 
of 852 spectral channels. Random Gaussian noise is then added 
to each of these channels at a level appropriate to the instrument 
sensitivity being simulated. 

At the same time as the data cube is generated, we also pro- 
duce a 'truth catalog' for the field which is comprised of 120//m 
fluxes (to allow comparison with the sensitivity of current and 
future instruments such as PACS and Spittzer), redshifts, posi- 
tions and source classifications for objects in the 128" x 128" 
field of view. This truth catalogue is also used to calculate the 
conventional '20 beams per source' confusion limit which, for 
this field, is 4.3mJy. 
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Fig. 1. Spectral energy distributions used for the models 
From the top downwards these SEDs correspond to following SED 
types in the models: Normal cold, Normal, starburst M82, starburst, 

10" L FIR LIG, 10" 5 L FIR LIG, Cold LIG, Hot ULIRG, Cold ULIRG, 
Syl, Sy2. 



3. Results of Models 

To determine the completeness of source identification via far- 
IR line emission, we search the simulated data cube for sources 
with strong line emission; we then compare the resulting list of 
sources with the 'truth catalog' derived from the original simula- 
tions. A line detection is conducted in the following way: at each 
image pixel the continuum is removed by fitting and subtracting 
a 4th order polynomial; the standard deviation, <x, of the result- 
ing continuum-subtracted spectrum is then calculated; a source 
is defined to contain a line if there exist channels in the spectrum 
in which the line flux is greater than some pre-defined multi- 
ple of this rms noise, ie. ncr. We find through comparison to the 
'truth catalog' that a reliable, ie. no false detections, list of line 
emitting sources is produced when n = 4.5, ie. when lines with 
significance > 4.5<x are selected. 

Figure 2 compares the broad band image of the field at 
120yum (left) with a continuum-subtracted image at a wavelength 
matching the faintest detected line-emitting source, at a wave- 
length of ~209yum (right), for a simulation with Icr noise of 
0.4mJy. We note that this line-emitting source has a continuum 
flux of 0.56mJy, nearly an order of magnitude fainter than the 
traditional confusion limit. However it is still significantly de- 
tected through its line emission. The number of sources recov- 
ered through their line emission as a function of continuum flux, 
compared to the initial distribution of fluxes in the input catalog 
is shown in Fig. 3. This clearly demonstrates the power of the 
technique, in that we recover all sources above the conventional 
confusion limit and a substantial fraction of those with fluxes up 
a factor of 10 below it. As expected, the fraction detected drops 
as the source gets fainter. The potential for breaking the con- 
tinuum confusion limit through the detection of sources in their 
far-IR line emission is clear. 

It is possible not only to locate a source in RA-DEC via 
blank-line surveys but also, very importantly, to determine its 
redshift. This requires several assumptions to be made. In some 
cases we detect more than one emission line in a given object, 
making redshift determination quite simple. In most cases, how- 





Fig. 2. Source detected in line below confusion limit 
Left: Continuum 120//m image of the field. Right: Narrow band image 
of the field at the wavelength at which a 560//Jy 120//m continuum flux 
source is detected by its line emission. Note that the source is clearly 
detected in the line image (the rightmost bright pixel, at location (12, 
8)) while in the continuum image nothing can be seen because of the 
confusion noise. Greyscales not the same, but chosen to make this effect 
clear. The additional source seen to the left in the line image is below 
our 4.5 <r detection threshold. 



ever, only one line is detected: in such cases the single lines de- 
tected are usually one of two bright far-IR emission lines [CII] 
157.7/mi or [OI] 63.2//m. In a typical dust SED which is charac- 
teristic of dust at a temperature of ~35K, the [CII] line lies long- 
ward of the SED peak, while [OI] lies at shorter wavelengths. If 
we assume that the dominant contribution to the confused con- 
tinuum emission comes from the line emitting source, then by 
considering the position of the line relative to the SED peak it 
is possible to guess which far-IR species is responsible for the 
emission line and thus get an accurate determination of redshift. 
The efficacy of this is shown in Figure 3, where the inferred red- 
shifts of detected sources with S m> m > lmJy are compared to 
the catalogue as a whole. As can be seen, we recover 100% of 
the redshifts for sources at z<2.5 that are brighter than lmJy 
ie. ~5 times fainter than the conventional confusion limit. By 
z~2.5, both of the two strongest FIR lines, [CII]157.7yum and 
[OI]63.2/im have moved out of the passband of our simulated 
instrument. We therefore do not expect to obtain any of these 
higher redshifts in the current simulation. 

We thus find that blind, blank-field spectral surveys are capa- 
ble not only of detecting sources well below the confusion limit 
but of obtaining their redshifts at the same time. Early analysis 
suggests we can distinguish between different luminosity evolu- 
tion models (the burst and bright models from Pearson (2006)) 
with a survey 10'xlO' or less in size using these techniques. 
This is something not possible using broadband number counts. 
Operation at higher redshifts will come from additional narrow 
lines not in the current model, such as the strong [SiII]34.8yum, 
[SIII]33.5/mi and [OIV] 24.9yum lines, some of which can be 
used as emission line diagnostics. Template correlation working 
on redshifted PAH features will also allow us to probe high red- 
shifts. 



4. Conclusions 

We have examined the potential for of imaging spectroscopic in- 
struments in studying the sources responsible for the CIB. The 
main limiting factor in CIB studies for far-IR missions is con- 
fusion resulting from the high continuum source density and 
large instrumental beams. The use of imaging spectrometers en- 
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Fig. 3. Recovery of Objects, with redshifts, from the FIR data 
cube 

Top: Shows the input flux distribution at 120//m from the 'truth catalog' 
(dotted line) and the sources detected from their line emission as a func- 
tion of their 120/mi continuum flux (solid line) for a 4.5<x line detection 
threshold and a lcr noise of 0.4mJy per channel. The conventional con- 
tinuum confusion limit of 20 beams per source is shown as a vertical 
line. Note that many sources below this confusion limit are detected, 
with sources having continuum fluxes nearly an order of magnitude be- 
low the limit being accessible. Bottom: The redshift distribution from 
the 'truth catalog' for sources brighter than lmjy at 120/jm (dotted line) 
and the redshift distribution for sources brighter than lmjy at 120/jm re- 
covered from the 4.5<x line detected sources. Note that all sources with 
z < 2.5 and F i2 o > lmjy have been recovered. 



ables the detection of sources by their line rather than continuum 
emission, allowing us to detect sources whose continuum flux 
is up to an order of magnitude below the conventional confu- 
sion limit. Furthermore, a by-product of the analysis of the line 
and continuum emission of sources detected in this way can be 
used to determine source redshifts. We believe there is consid- 
erable mileage in deep, blank-field spectral line surveys in the 
far-infrared and will discuss the results of more extensive mod- 
eling in a future, more substantial paper. 
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